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CPT invariance is one of the most fundamental theo- 
retical concepts; its violation would have a serious impact 
on physics in general, and would require new physics be- 
yond the standard model (SM). CPT violation requires 
the breakdown of some fundamental underlying physical 
assumption in the new physics beyond the SM, for exam- 
ple, violation of Lorentz invariance [l| . Several searches 
for CPT violation have been carried out; for example, 
the Belle and BaBar collaborations have published mea- 
surements of CPT-violating parameters in the neutral B- 
meson system and the CPLEAR, KLOE, and KTeV 
collaborations have done so in the neutral K-meson sys- 
tem 013. 



In the presence of CPT violation, the flavor and mass 
eigenstates of the neutral B mesons are related by \Bl) ~ 
p VT^ \B°) + qVm\W) and \Bh) = p^/T + I|pO) - 
gVl ~ z\B^), where \Bl) {\Bh)) is a light (heavy) mass 
eigenstate. Here z is a complex parameter accounting for 
CPT violation; CPT is violated if z 7^ 0. In the decay 
chain T{4S) B°B° /rcc/tag, where one of the B- 
mesons decays at time tree to a reconstructed final state 
/rcc and the other decays at time ttag to a final state 
/tag that distinguishes between B^ and the general 
time-dependent decay rate with CPT violation allowed 
is given by [1] 
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whcrc Vd = (r^f 

rriL, At = trec-ttag, and the •^Bo,BO-^f,^^,f^, 
evant decay amplitudes. If free is a CP eigenstate {fcp), 
a parameter Xcp, which characterizes CP violation, can 
be defined as Xcp = {qlv){-^B«^fcpl ^B^^fcp)- The 
SM predicts |Acp| — 1 and Im{ric p\c p) — sin2(/)i for 
the case fcp = J/i^K^, where ?]cp is the CP eigenvalue 
of the final state. 

In this paper we report improved results on the CPT- 
violating parameter z and on the normalized total-decay- 
width difference ATd/Td in B° J/ipK^ {K° = 
K°,Kl), D'^*^-h+ {h+ = 7r+ for D" and 7r+,p+ 

for £)*-), and 5° ^ D*-i+vt {£+ = e+,fi+) decays g. 
Most of the sensitivity to TZe{z) and ATd/Td is obtained 
from neutral B- meson decays to fcp- while Im(z) is 
constrained primarily from other neutral i3-meson decay 
modes. 

The data sample of 535 x 10® BB pairs used in 
this analysis was collected with the Belle detector at 
the KEKB asymmetric-energy e+e" collider Q (3.5 on 
8.0 GcV) operating at the T(45') resonance. The T{4S) 
is produced with a Lorentz boost of p-f = 0.425 along 
the Z axis, which is antiparallel to the beam direc- 
tion. Since BB pairs are produced approximately at rest 
in the T(4iS') center-of-mass system (cms), At can be ap- 
proximated from AZ, the difference between the Z coor- 
dinates of the two B decay vertices: At ~ AZ / {/3jc). 

The Belle detector is a large-solid-angle magnetic 
spectrometer that consists of a silicon vertex detector 
(SVD), a 50-layer central drift chamber, an array of 
aerogel Cherenkov counters, a barrel-like arrangement of 
time-of-flight scintillation counters, an electromagnetic 
calorimeter comprised of CsI(Tl) crystals, located inside 
a superconducting solenoid coil that provides a 1.5 T 
magnetic field. An iron flux-return located outside of the 
coil is instrumented to detect mesons and to identify 
muons. Two inner detector configurations are used; a 
2.0 cm radius beam pipe and a 3-layer SVD are used for 
the first data set (DS-I) of 152 x lO'' BB pairs, while a 
1.5 cm radius beam pipe, a 4-layer SVD, and a small- 
cell inner drift chamber are used to record the remaining 
data set (DS-II) of 383 x 10® BB pairs. 

We reconstruct — > f^cc decays in the 
D-Tr+, D*-TT+, D*-p+, and D*^i+ve channels. We also 
reconstruct J/ipK^ and D^tt"*" to precisely deter- 



mine parameters for the At-resolution function model in 
neutral B decays. For the J/ipKg and J/tpKj^ modes, 
we use the same selection criteria as in Ref. |ll| . Candi- 
date J/ipK~^ events are selected from combinations of a 
charged track and a J/ip candidate using the same selec- 
tion criteria as in J/ipKg. Charged and neutral charmed 
mesons are reconstructed in the — > K^tt^tt^ and 
— >■ K~^Tr~ , K'^n~TT'^ , K~^Tr~TT~^TT~ decay modes, re- 
spectively. The invariant mass of their daughters, Mkutt, 
is required to be within 45 MeV/c^ (^5(t) of the nominal 
ZJ-meson mass for the mode with 7r°, or 30 McV/c^ (^6(t) 
for the other modes. Candidate D*~ mesons are recon- 
structed in D^t:~ combinations, in which the mass dif- 
ference Mdiff between the D*^ and D° candidates is re- 
quired to be within 5 MeV/c^ (^8ct) of the nominal value. 
Candidate p+ mesons are reconstructed from tt+tt® com- 
binations with invariant mass within 225 MeV/c^ of the 
nominal p"*" mass. The D*~ candidates for the final state 
D*~£~^vg are reconstructed using the D*~ and D° decay 
modes listed above, where the detailed selection criteria 
arc described in Ref. (T2j . 

We identify 5° or 5+ candidates in modes other than 
B'^ — > J/ipK^ or D*^£^i^g using the beam-energy con- 

p, and the en- 



strained mass, Mbc = \/ {E^ 
ergy difference, AE = E^ 



beam 

^beam> ^^cre E^^^^^ is the 
beam energy in the cms, and E^ and p*g are the cms 
energy and momentum of the reconstructed B candi- 
date, respectively. The signal region for the Mbc is de- 
fined as 5.27 GeV/c2 < M^c < 5.29 GeV/c^ for aU de- 
cay modes, while that for AE is decay-mode dependent: 
I AS I < 40 MeV for J/^pK^ and J/'ipK+; |A£;| < 45 MeV 
for 0-7:+; \AE\ < 70 McV for D*~n+- -50 MeV < 
AE < +80 MeV for D*' p+ , and |A£;| < 60 McV for 
Z?"7r+. Candidate 5° J/ipK^ decays are selected 
by requiring 0.20 GeV/c < < 0.45 GcV/c. For 

B^ — )■ D*~t^vi decays, the energies and momenta of the 
B meson and D*l system in the cms satisfy M^^ = [E'^ — 
Eo'-i^? - iW + \Pp'i? - 2|p-||p-.,|cos0B,D.,), 
where M,^,, is the neutrino mass and cosOb.d-i is the 
angle between ■p*g and p^*^. We calculate cosOb,d*i set- 



ting Afyj = and E*^ = E^^^^. The signal region is 
defined as |cos6'b^d<.^| < 1.1. In the cos9B,D*t signal 
region, B^ — >■ D**~£^vi decays are also reconstructed. 
Since the At distribution is expected to be the same as 
that in 1/^, we treat B^ D**~g~^vi decays as 



4 



signal. 

The event-by-event signal and background probabili- 
ties are estimated from signal and background distribu- 
tions of the kinematic parameters, Aibc: AiJ, \p^\, and 
cos9B,D'i- Signal and combinatorial background distri- 
butions in A/bc are modeled by Gaussians and an em- 
pirically determined background shape with a kinematic 
threshold originally introduced by ARGUS [H], respec- 
tively, while those in AE arc modeled by the sum of two 
Gaussians and a first-order polynomial, respectively. The 
model parameters for the signal and combinatorial back- 
ground distributions in the J/ipK^ and J/ipK~^ modes 
are determined from a two-dimensional fit to the Mbc- 
AE distributions in data. In Monte Carlo (MC) simu- 
lations of the /i+ and £)°7r+ modes, in addition to 
combinatorial background, we find a background contri- 
bution, which comes from charged and neutral i3-mcson 
decays with one or more particles missed in their recon- 
struction, and which peaks in Mbc (peaking background). 
The model parameters for the signal and combinatorial 
background distributions are determined from signal and 
sideband A/bc distributions in data, while those of the 
peaking background are modeled by an ad-hoc distri- 
bution obtained from MC simulation. In addition to 
the combinatorial background, we find from MC simu- 
lation that the background in J /ipK^ is mainly com- 
prised of {cc)K^ modes except for contributions from 
J/tl'Kl, J/ipK^TT^, J/-07r°, and charged i?-meson de- 
cays. For D*^£'^i'i, there is an additional background 
from D** i^iyg. For the J/tpKl and D*~(.^Vi modes, the 
signal and noncombinatorial background distributions in 
and cos0B,D'-e, respectively, are modeled using MC 
simulation, while the combinatorial background distribu- 
tions are obtained from sideband regions of the J/ip and 
D*~ , respectively. 

The 6-flavor of /tag is identified from inclusive proper- 
ties of particles that are not associated with the ^■ 
/roc decay. The tagging information is represented by two 
event-by-event parameters, the 6-flavor charge gtag and 
an MC-determined flavor-tagging dilution factor r p^ . 
The parameter r ranges from r = for no flavor discrim- 
ination to 7' = 1 for unambiguous flavor assignment. For 
events with r > 0.1, the wrong tag fractions for six r 
intervals, wi (I =■ 1 . . .6), and their differences between 
B° and B° decays, Awi, are determined using the data 
sample as described later. If r < 0.1, we set the wrong 
tag fraction to 0.5 so that the event is not used on flavor 
tagging. 

The vertex position is reconstructed using charged 
tracks that have sufficient SVD hits |l^ . The /roc vertex 
for the modes with a J /tp is reconstructed using lepton 
tracks from the J /if) decay, while in modes without a 
J 1^ the /rec vertex is reconstructed by combining the 
or D~-meson trajectory and the remaining charged 
track forming the B-meson candidate; the slow 7r~ from 
the D*~ decay is not included because of its poor po- 
sition resolution. The /tag vertex is obtained from se- 
lected well-reconstructed tracks that are not assigned to 



TABLE I: Number of events A'^ov and purity in the signal 
region for each decay mode. 



D decay mode 


AT 


Jr^Ul ity {/oj 


■J iW^s 


771 


Q7 n 




10966 


59.2 


D-TT + 


39366 


83.2 


D'-Tr+ 


46292 


81.5 


D*-p+ 


45913 


66.3 




383818 


75.2 


J/ipK+ 


32150 


97.3 


D°n+ 


216605 


63.9 



/roc- A constraint on the interaction region profile (IP) 
in the plane perpendicular to the Z axis is also applied 
to both /roc and /tag reconstructed vertices. We model 
the resolution function R{At) as a convolution of four 
sub-components p^ : detector resolutions for f^cc and 
/tag vertex reconstruction, boost effect due to nonpri- 
mary particle decays in /tag, and dilution by the kine- 
matic approximation At ~ AZ /{/S^c). Nearly all model 
parameters are determined using the data as described 
later. The exceptions are the parameters for the boost 
effect and kinematic approximation, which are obtained 
using MC simulation. For candidate events in which both 
B vertices are found, for further analysis, we only use 
events with vertices that satisfy ^roc < 250, ^tag < 250, 
and \At\ < 70 ps, where ^roc (^tag) is the of the /roc 
(/tag) vertex fit calculated only along the Z direction [13 |. 

After flavor tagging and vertex reconstruction, we 
count the number of events remaining in the signal re- 
gion A"ov and estimate the purity for each decay mode. 
The values of TVov and purity for each mode are listed in 
Table H 

We determine three major physics parameters TZe{z), 
Tm{z), and ATd/Td together with five other physics 
parameters Tgo , Tg+ (neutral and charged B-mcson 
lifetimes), Anid, |Acp|, and a.vg{ricp^cp) in a simul- 
taneous 72-parameter fit to the observed At distri- 
bution. The remaining 64 parameters are the At- 
resolution function model parameters (34), flavor-tagging 
parameters wi and Awi (24), and background param- 
eters for _B" D*~t^vt (6). The nonphysics pa- 
rameters arc determined separately for DS-I and DS- 
II. An unbinned flt is performed by maximizing a like- 
lihood function defined by L{TZe{z),Tm{z), AVd/Td) ~ 
n, i*(7^e(^), Jm(z), AFd/Fd; Af , gj^J, where the prod- 
uct is over all events in the signal region. The likelihood 
for the i-th event U is given by 

L^ = {1- .foi).fL,V{Af-, fl,, flj ® R^{Af) 

+ (1 - /oi) E /b'i;^bkg(Af ) + /oiPoi(Af ). (4) 

k 

The first term accounts for the signal component, where 
/gi is an event-by-event signal fraction. In Eq. (j4]) V 
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is modified from Eq. ([T]) by including tlie event-by-event 
incorrect-tagging effect, wl and AwJ, and the symbol ® 
indicates a convolution with the At-resolution function 
W{At). The second term accounts for the background 
component, where f^^^ is an event- by- event background 
fraction and k runs over all background components. 
The signal and background fractions arc normalized to 
/sig + Sfc /bkg = 1- The At-distribution for the combina- 
torial background component is modeled using the side- 
band region of Ai?-Mbc, \Pb\, or cos6b,D'1 space, while 
the At distribution for the peaking-background compo- 
nents are modeled by MC simulation. The third term ac- 
counts for a small but broad (At outlier) component that 
cannot be described by the first and second terms, where 
/oi is an event-dependent outlier fraction and Pol (At) 
is a broad Gaussian. In the nominal fit, we account 
for Cabibbo-Kobayashi-Maskawa-favored S — > Z? decay 
via h — > cud (CFD) but neglect the contribution from 
Cabibbo-Kobayashi-Maskawa-suppressed B ^ D decay 
via h — ucd (CSD) both in /roc and /tag- The effect of 
the CSD is included in the systematic uncertainty. 

From the fit to the data, we obtain TZe{z) = (4-1.9 ± 
3.7) X 10^2^ Im{z) = (-5.7±3.3) x IQ-^, and AVd/Td = 
(-1.7±1.8)xl0"^ together with Tso = 1.531±0.004 ps, 
tb+ = 1.640 ± 0.006 ps, Ama = 0.506 ± 0.003 ps^^, 
|Acp| - 1 = (1.1 ± 3.8) X 10-3, and arg{vcpXcp) = 
—0.700 ± 0.042, where all uncertainties are statistical 
only. The fit has a twofold ambiguity in the sign of 
'R-e{r]cpXcp), 'Tie{z) and AVd/Td change signs depend- 
ing on its sign. We take the solution with positive 
T^e{ricp^cp), which is the result of the global fit fl6l |. 
The correlation coefficients p between two of the three 
major physics parameters are PTie{z) ,Xm{z) = —0.17, 
P7^e(^),A^d/^d = +0.08, and Pi,n{z).Ard/ra = +0-09. The 
largest correlation coefficient between a major physics 
parameter and any other fit parameter is P'R,e{z),Amd ~ 
-1-0.24. The fitted values of |Acp| and arg(?7cpAcp) give 
sin 201 = 0.645 ± 0.032(stat), which is consistent with 
our dedicated sin 2(f)i measurement with the same data 
sample 
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because the major physics parameters are 
consistent with zero. Figures [1] and [2] show the At distri- 
butions for events with /rec = J/ ipK'^ cases and the other 
cases, respectively, with the fitted curves superimposed. 



To illustrate the CPT sensitivity of our measurements, 
we plot the deviations of the asymmetries from a refer- 
ence asymmetry obtained from the nominal fit parame- 
ters but setting 7^e(z) = Im{z) = ATa/Td = in Fig. El 
where (a), (b), and (c) show those for CP asymmetries 
of — > J/ipKg, J/i/jK1, and opposite-flavor i3-meson 
pairs, respectively; (d) shows asymmetries between the 
opposite-flavor and same-flavor i?-mcson pairs. Asym- 
metries are obtained from events in At bins without 
background subtraction, where the events are required 
to have r > 0.5. We superimpose the deviations of the 
asymmetries for the nominal fit curves and those with 
one parameter shifted by ~ 5 times the statistical uncer- 
tainty in each subsample fit. For illustration, the most 
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FIG. 1: At distributions for events with flavor tag quality 
r > 0.5, where (a) and (b) correspond to free = J/4>K% and 
J/rpKl cases, respectively. Events are separated according to 
tagged /tag flavor, where the solid and dashed curves are for 
9tag = -1-1 and —1 events, respectively. The two chain curves 
below the flt curves indicate the sum of the background and 
At-outher components for each /tag flavor, which are almost 
indistinguishable because of their similar shapes. 
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FIG. 2: At distributions for events with flavor tag quality 
r > 0.5, where (a, b) and (c, d) correspond to flavor-specific 
/,g^ = B° and W cases, respectively. The dashed curve below 
the solid fit curve is the sum of the background and At-outlier 
components. 



appropriate parameter is chosen in each plot. 

Table |ll] lists the systematic uncertainties on the ma- 
jor physics parameters. The total systematic uncertainty 
is obtained by adding the contributions in Table |TT] in 
quadrature. The dominant contributions arc from the 
tag-side interference (TSI) [iTj and vertex reconstruction; 
the next largest contributions are from fit bias. 

The TSI effect arises from the interference between 
CFD and CSD amplitudes in /tag. In general, the pres- 
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(a) so — J/yKg decay 



(b) SO — J/yKp decay 
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(c) Flavor-specific S decays 
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(d) Flavor-specific S decays 
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FIG. 3: Deviations of the asymmetries from tlie reference 
asymmetry. The crosses with error bars are data. The solid 
curves are deviations for the nominal fits. The dashed curves 
are with 7^e(^) = +0.28 for (a) and (b), Im{z) = -0.03 for 
(c), and AFd/rd = —0.16 for (d) (see text for details). 



TABLE II: Summary of systematic uncertainties on the major 
physics parameters. 



Source 


S{TZe{z)) 


S{Im{z)) 




Vertex reconstruction 


0.008 


0.0028 


0.009 


At-resolution function 


0.003 


0.0004 


0.002 


Tag-side interference 


0.028 


0.0006 


0.001 


CSD effect 


0.004 


0.0008 


0.003 


Fit bias 


0.012 


0.0013 


0.005 


Signal fraction 


0.004 


0.0002 


0.002 


Background At shape 


0.005 


0.0001 


0.002 


Others 


0.001 


< 0.0001 


0.002 


Total 


0.033 


0.0033 


0.011 



ence of CSD introduces new terms in Eqs. ([2|) and ([3]) 



■^B»^f CSD 



%ecexp[i(+03 + <5/,,J] 
Rf^^^ exp [i (-03 + Sf^^J] 



(5) 



where and Sf^^^ are the mode-dependent ratio of the 
CSD to CFD amplitudes and the relative strong-phase 
difference between the CSD and CFD amplitudes, respec- 
tively, and 03 = 67.2° [l^. For the tag-side parameter, 
-^/tag 9,nd Sf^_^^ are "effective" values because /tag is an 
admixture of several decay modes, some of which do not 
have a corresponding CSD. The effective R/^^^ and (5/^^^ 



parameters are estimated using the D*~t^vi sam- 

ple [ill- We perform fits to the major physics parameters 
varying the terms from Eqs. ([SJ into Eqs. ([U and 
The deviation from the nominal fit is quoted as a sys- 
tematic uncertainty. 

The CSD effects in /i-cc are investigated by performing 
fits of the major physics parameters varying the and 
5f\^^ parameters introduced in Eqs. For the D~tt'^ 
and D*^n^ modes, we use Rdtt = 0.02 or Rd'-k = 0.02 
predicted in Ref. [l8|, and 5ui')h computed from mea- 
surements of CP- violating parameters in the relevant B 
decays [l^. We quote fitted deviations as the system- 
atic uncertainties. For the D*~ p'^ mode, we assume 
Rd'p = 0.02, and allow So-p to be 0°, 90°, 180°, or 
270°, because of the absence of CP- violating parameter 
measurements. Wc quote the largest fitted deviation as 
the systematic uncertainty. 

The systematic uncertainty due to vertex reconstruc- 
tion is estimated as follows. We repeat fits by changing 
various requirements or parameters used in the vertex 
reconstruction: the IP constraint, the track selection cri- 
teria, and the calibration of the track position and mo- 
mentum uncertainties. The deviation from the nominal 
fit is quoted as the systematic uncertainty. Systematic er- 
rors due to imperfect SVD alignment arc estimated from 
MC samples with artificially varied alignment constants. 
Effects from small biases in the AZ measurement ob- 
served in e^e^ — > /i^/i^ and other control samples are 
accounted for by applying a special correction function 
and including the variation from the nominal result into 
the systematic uncertainty. 

We estimate the fit biases (55,'''^ s , 5^^'^'', . , and 5^}^^ 

Ke[z)^ l.m[z)^ Aid/id 

using an analysis procedure with fully simulated MC 
samples. We generate sets of At distributions with statis- 
tics similar to data, fixing {TZe{z),Xm{z), ATd/Td) = 
(0, 0, 0) or varying one of the three input parameters to 
7^e(z) = ±0.01, Im{z) = ±0.01, or ATa/Td = ±0.05. 
We perform a full-parameter fit to each generated dis- 
tribution without the background component, and take 
deviations of the fitted three parameters from the input 
value as the bias. We quote the average value of biases 
in the above seven samples. These effects are included 
into the systematic uncertainty after symmetrization. 

The systematic uncertainty due to the At-resolution 
function is estimated by varying by ±2a each resolution- 
function parameter determined from MC, and repeating 
the fit to add each variation in quadrature. We also 
take the systematic effect from the At-outlicr elimina- 
tion criteria into account in the systematic uncertainty 
by varying each criterion and adding each variation in 
quadrature. 

The most precise previous results on the CPT- 
violating parameter and ATd/^d in the neutral B- 
meson system were obtained by the BaBar collabora- 
tion. They found 7ee(Acp/|Acp|)7ee(z) = +0.014 ± 
0.035(stat) ± 0.034(syst), Jto(z) = (-13.9 ± 7.3(stat) ± 
3.3(syst)) X 10-3, and sgn(7^e(Acp))A^d/^<J = -0.008± 
0.037(stat)±0.018(syst) [sB. For 7^e(Acp/|Acp|)7^e(z), 



7 



our result is (+1.5 ± 3.8) x 10~^, where the total error 
is quoted. Our result is consistent with Rcf. Q and im- 
proves the overall precision by factors of 1.3 to 2.0 for all 
parameters. 

In summary, we report a new search for CPT violation 
with an improved measurement of the CPT- violating pa- 
rameter z and normalized decay-rate difference AFd/Fd 
in ^ J/V-A'g, J/i^Kl, D~TT+, D*-n+, D*' p+ , and 
D*^l^vi decays using 535 x 10^ BB pairs collected at 
the T{AS) resonance with the Belle detector. We find 

ne{z) = [+1.9±3.7(stat)±3.3(syst)] X 10"^ 
Im{z) = [-5.7±3.3(stat)±3.3(syst)] X 10"^ and 
AVd/Td = [-1.7± 1.8(stat) ± l.l(syst)] X 10"^ 

all of which are consistent with zero. This is the most 



precise measurement of CPT-violating parameters in the 
neutral _B-meson system to date. 
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